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Abstract: We demonstrate the effectiveness of averaging over the Wilson parameter r 
(which has been proposed earlier) in removing the cutoff effects of naive Wilson fermions in 
both the anomaly term and the pseudoscalar density term in the flavor singlet axial Ward 
Takahashi identity at 0(g 2 ) involving slowly varying background gauge fields. We show 
that it is the physical fermion contribution which is largely influenced by the r averaging. 
We have studied the possible interplay between finite size and cutoff effects by investigating 
in detail naive, 0(a) improved and OStm Wilson fermion cases for a range of volumes and 
lattice fermion mass (am). For naive Wilson fermions r averaging is shown to remove 
the effects of the interplay. We have shown that for the pseudoscalar density term to 
0(g 2 ) the lattice result differs from the continuum result by exhibiting considerable am 
dependence which appears to be a manifestation of cutoff effects with naive Wilson fermion. 
The pseudoscalar density term to 0(g 2 ) is shown to be almost independent of am when 
r-averaging is performed. 
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1 Introduction 

Recently we have studied [1] the emergence of the chiral anomaly in the continuum chiral 
limit with Osterwalder-Seiler twisted mass (OStm) Wilson fermions [2] in lattice QCD 
in comparison with naive [3], 0(a), and 0(a 2 ) improved [4-6] Wilson fermions. Here a 
denotes the lattice spacing. Karsten and Smit [7], Kerler [8, 9] and Seiler and Stamatescu 
[10] have demonstrated some time ago that naive Wilson fermions on the lattice can indeed 
reproduce the chiral anomaly in the infinite volume chiral limit in the continuum. However 
for a lattice simulation, understanding the effects of the finite volume, finite cutoff and 
nonzero fermion mass is very crucial. Possible interplay among various lattice artifacts in 
the unimproved theory and how various suggested improvements help to ameliorate the 
situation also needs to be investigated. Lattice calculation of the chiral anomaly to one loop 
perturbation theory with finite lattice fermion mass (am) and finite box size (L) provides 
an excellent laboratory to address some of these issues in a simple setting. 

Aoki [11] proposed taking the average over the Wilson parameter r = ±1 for improving 
the scaling behavior of chiral condensate. David and Hamber [12] also proposed similar 
ideas. For a detailed discussion of the symmetries associated with r — > —r and their 
implications, specifically that the r averaged correlators are affected by only 0((a 2 ) k ), k = 
1,2,..., discretization effects, see Frezzotti and Rossi [13]. Twisted mass (tm) fermion 
[14-16] has an effective r averaging, but at finite lattice spacing there are isospin breaking 
cut off artifacts. However, in all simulations carried out until now, 0(a 2 ) effects have been 
shown to be under control and to disappear in the continuum limit. Large 0(a 2 ) effects 
are only visible in the neutral pion mass [17]. This observation has been discussed in detail 
and explained theoretically in Ref. [18]. 
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We note in passing that since positivity property of the quark determinant is not 
satisfied by the OStm Wilson (W) QCD action, numerical simulation is only possible with 
a mixed action with OStm WQCD action for the valence quarks and tm WQCD action for 
the sea quarks. 

r averaging of lattice QCD with naive Wilson fermions looks interesting. The flavour 
singlet axial anomaly in this case can be studied analytically under certain assumptions and 
the effect of the r averaging investigated quantitatively for each of the terms on the right 
hand side of the Ward Takahashi identity. In this work there is overwhelming evidence for 
improvement with respect to cutoff effects, once r averaging is done. In fact the r averaged 
result is much better than 0(a 2 ) improved Wilson fermion with tree level coefficients [4]. 
In addition there are interesting issues of interplay between 0{a) cutoff effects and finite 
volume effects, and its removal after r averaging. Contributions from the physical fermion 
and the different doubler fermions are also considered separately. Finally the effect of r 
averaging on the pseudoscalar density (PSD) term on the right hand side of the Ward 
Takahashi identity is also studied. In the following paragraphs of the introduction we 
briefly discuss some of these issues. 

Since the emergence of anomaly is intimately tied with the removal of fermion doublers, 
the behavior of their contributions as one varies the parameters of the theory is of interest. 
Recall that as the Wilson parameter r is set to zero, the anomaly contributions from the 16 
species (the physical fermion and the 15 doublers) all cancel each other. Karsten and Smit 
showed that for non-vanishing r, anomaly is independent of r in the limit a — s> 0. However, 
on a lattice with finite cutoff (practical situation), doubler masses are not exactly infinitely 
large and they do contribute to the anomaly depending upon the value of r. We study the 
doubler contribution [4, 7] as a function of r and am. 

We also compare the finite volume dependence of the one loop calculation of the 
anomaly with naive Wilson action and the improved actions. The effect of small momentum 
behavior of the integrand as the fermion mass is lowered is investigated with both Periodic 
(P) and Anti-Periodic (AP) boundary conditions (BC). We find that the naive Wilson 
fermion has larger finite volume effects than the improved ones possibly because of the 
interplay between the finite volume and cutoff effects. Even though OStm Wilson fermion 
has a better chiral behavior as far as approach to the chiral limit is concerned than the 
0{a) improved Wilson fermion in the infinite volume limit [1], we find that with PBC, 
finite volume effects are more pronounced with the former than with the latter. APBC 
improves the convergence behavior as L becomes large for both naive and OStm Wilson 
fermions. The most striking result, however, is that finite volume effects of naive Wilson 
fermions are largely removed with r averaging, indicating that finite volume effects of naive 
Wilson fermions observed in this calculations are mostly at 0(a). 

At the tree level, in the continuum limit, the flavour singlet axial vector current is con- 
served in the chiral limit because both the PSD term and the contribution from the Wilson 
term (denoted by x later) vanish. To 0(g 2 ), in the continuum limit, it can be seen that the 
PSD term is independent of the fermion mass for a slowly varying background gauge field 
[19] and produces the negative of the anomaly term, leading to a cancellation. However 
note that this cancellation is due to the assumption of slowly varying background gauge 
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field and at this order only. For general gauge fields and to any order of g 2 , cancellation 
does not occur and the flavour singlet axial vector Ward Takahashi identity is anomalous. 
For 0(a) improved Wilson fermions, at 0(g 2 ) the cancellation has been demonstrated in 
the chiral limit for constant background gauge field by Hamber and Wu [4]. Our starting 
point is the flavour singlet axial vector Ward Takahashi identity on the lattice with naive 
and improved Wilson fermions. Our emphasis is to study the effect of r averaging on the 
anomaly and the PSD terms individually, as discussed above. In the case of naive Wilson 
fermions, both the PSD and the anomaly terms individually show strong lattice mass de- 
pendence. By using the Reisz's power counting theorem [20] for the PSD term to 0(g 2 ) 
for a slowly varying background gauge field we show that the PSD term is independent of 
fermion mass in the continuum limit. This result indicates that the observed mass depen- 
dence with naive Wilson fermion is a pure lattice artifact which is largely removed by r 
averaging leading to remarkable am-independence. 



2 Naive, 0(a) improved and OStm Wilson Fermions 

The flavor singlet axial Ward Takahashi identity for the naive Wilson fermions on the 
Euclidean lattice reads 



(A b J 5fl (x)) = 2m(i{j xl M + ( X a 



(2.1) 



where (O) denotes the functional average of O for a background gauge field. Explanation 
of other terms are as follows: 



The backward derivative, A^f(x) = - [f(x) — f(x — fi)] , 

and ( Xx ) = -Tr ace[~/ 5 {GW + WG)} . 



(2.2) 



The Green function 



G(x,y) = (x 



[^D^ + W + m] 



y) , 



(2.3) 



where 



'Pixy 2a 



_ Tn A 



u 



(2.4) 



For ease of comparisons and to facilitate discussions, here we collect the expressions 
for chiral anomaly in the naive, 0{a) improved and OStm Wilson fermions [1]. 
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The contribution to axial vector Ward Takahashi identity from the Wilson term in the 
unimproved theory 



(Xx) 



2 g 2 e^px trace F^(x)F pX {x) -^-^ ^ cos(p /J a)cos(p^a)cos(p 



x W (p)\cos(p x a)[m + W (p)\ - 4 - sm 2 ( Px a) (Goip))' 6 , 
l a 

g 2 

~~ VuXi e ^P x trace Ffiiy{x)F p x(x)I(am,r, L) 



(2.5) 
(2.6) 



where 



W {p) = -^[(l-cos(a^))] , 

/ y r ^ 

Goip) = [—^2 sil)2 ( a Pn) + ( m + _ ~ cos ( a ^)]) 2 



(2.7) 



and trace is in color space. 

Explicitly, Y, P = (r) 4 Eni,n 2 ,n 3 ,n4 where ™i, n 2 , n 3 , n 4 = 0, 1, 2, 3, • • • . In the infinite 
volume, chiral limit, / — > — ^- so that —\l = I — > 1. In all our plots it is the function 
I(am, r, L) called anomaly integral which we have plotted. 

For 0{a) improved Wilson fermion 



(xl) 



2 g 2 e^px trace Fp V (x)F pX (x) J^Xi ^ cos(ap M ) cos(ap p ) cos(ap A ) [Gq 



cos(a Pu )[m + W (p) + W£(p)] 

r 1 i 

- 4 - sin(ap u )(sm(ap u ) - -sm(2ap u )) [W (p) + Wq(p)] . 
a £ - 



(2.8) 



Here 



Woip) + Wtf (p) = [J2 [^(1 - cos(app)) + ^(-1 + cos(2ap M )) 
g o(p) = f 4 E sin2 ( a ^) + ( m + E& 1 " cos («^)] + it" 1 + cos^)]]) 2 ^ (2.9) 



In the case of the OStm Wilson action 



(X° x s ) = 2 g 2 €p Up x trace Fp U (x)F pX (x) — ^ ^ cos(^a)cos(p !/ a)cos(p p a) 



x W (p) co S (p x a)W (p) - 4 - sin^(p A a)|(^ s (p)) 3 



(2.10) 



where 



(G 



(2.11) 



- 4 - 




Figure 1. Anomaly integral for 0.001 < am < 0.1 for r averaged naive Wilson fermions at L = 60 
for r = ±1 compared with OStm, 0(a) and O{o?) improved Wilson fermions with r = 1. 

3 Wilson parameter r Averaging 

In Fig. 1 we compare the anomaly integral for r averaged naive Wilson fermion with the 
OStm, 0(a) and 0(a 2 ) improved Wilson fermions. The r averaged anomaly integral for 
naive Wilson fermion shows better chiral behaviour compared with OStm Wilson fermion 
for am > 0.01. The cutoff effects are known to be proportional to even powers of a for both 
r-averaged and OStm Wilson fermions but they are not the same as illustrated in Fig. 1. 
Although the leading cutoff effect in r averaged naive Wilson fermion and OStm Wilson 
fermion is of 0(a 2 ), they have much better chiral behaviour than the 0(a 2 ) improved [4] 
Wilson fermion which employs tree level coefficients. For further improvement in the latter 
case one needs to tune the coefficients while OStm and r-averaged Wilson fermion do not 
need any tuning. 

In the following sections we present results of r averaging for the anomaly integral for 
different lattice volumes and am and contributions to the integral from different doubler 
sectors. 

4 Anomaly: contribution from different regions of the Brillouin zone 

It is instructive to separately consider the contribution to the anomaly integral from dif- 
ferent regions of the Brillouin zone in order to compare and contrast the contribution of 
the physical fermion and the doublers as a function of the different parameters. This has 
been studied in the chiral limit by Hamber and Wu [4]. However, for practical simulations 
the behavior at non-zero quark masses and nonzero lattice spacing is more relevant and we 
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study this issue in this section. Following Karsten and Smit, the limits on the momentum 
sum are changed from (— tt/cl, +ir/a) to (— 7r/(2a), 37r/(2a)) and further the momentum 
sum hypercube is divided into 16 smaller hypercubes corresponding to (— 7r/(2a), +7r/(2a)) 
and (+7r/(2a), +37r/(2a)) for each pu,/i = 1,2,3,4. Thus the total anomaly contribution 
is decomposed into the contributions from the five species-types and the Anomaly Inte- 
gral I = Iq — 4ii + 6/2 — 4/3 + I4. In Figs. 2(a) and 2(b) we compare the contributions 
in the region 0.01 < am < 1.0 for r = 1 for naive Wilson fermions and OStm Wilson 
fermions respectively. As expected, as the physical fermion gets lighter, its contribution to 
the anomaly is dominant compared to the doublers. Nevertheless, doublers' contribution 
is not insignificant. The first doubler contributes the most since it has the lowest mass 
among the doublers. Also note that the lowest mass doubler has more am dependence in 
the case of naive Wilson fermion compared to OStm fermion. 

To verify that the doublers do decouple when their masses become very heavy so 
that the sole contribution to the anomaly comes from the physical fermion, we study the 
behavior of Ii,i = 1,2,3,4 as the Wilson parameter r is increased. We realize that the 
reflection positivity is not guaranteed when | r \> 1, but in the present context we adjust 
r only to make the doubler masses heavy. 

In Fig. 2(c) we present the total contribution of the doublers to the anomaly integral for 
0.01 < am < 0.1 for naive Wilson fermions at L = 40 for 1 < r < 2. As expected, the total 
doubler contribution diminishes as r increases and it is not very sensitive to am in the range 
0.01 < am < 0.1. Fig. 2(d) shows, for r = ±1, total doubler contribution and physical 
fermion contribution separately. For each r, total contributions from the doublers have 
slight am dependence. Although for each r, the contribution from the physical fermion 
depends significantly on am, the figure clearly indicates that, on r averaging, the am- 
dependence will cancel out. 

5 Volume Dependence 

Since numerical simulations are performed at finite volume and finite lattice spacing, it 
is of interest to study the finite volume dependence of the anomaly integral as a function 
of the quark mass. Since finite size artifacts are different for the three fermion actions 
considered in this work, and finite volume and finite cutoff (1/a) effects may interfere with 
each other, the anomaly integral provides us an opportunity to explore this issue in this 
simple setting. 

First we discuss the case of PBC in all four directions. In Fig. 3(a) volume dependence 
of the anomaly integral for 0.01 < am < 0.1 for naive Wilson fermion for r = 1 and PBC 
is presented. We note that the behavior qualitatively changes as am changes from large to 
small. Moreover, irrespective of the quark mass, approximate L independence is achieved 
only for L > 40. In Fig. 3(b) the corresponding behavior for 0(a) improved Wilson fermion 
is shown. Here the L dependence is qualitatively the same across the range of am, and 
furthermore, L independence is achieved for L > 20. Volume dependence of the anomaly 
integral for OStm Wilson fermion is presented in Fig. 3(c). Here again the L dependence 
is qualitatively the same across the mass range but L independence is achieved only for 
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Figure 2. (a) Contributions from different regions of the Brillouin zone to the anomaly integral 
for r = 1 and 0.01 < am < 1.0 for naive Wilson fermions at L = 40 compared with the total 
contribution, (b) Contributions from different regions of the Brillouin zone to the anomaly integral 
for r = 1 and 0.01 < am < 1.0 for OStm Wilson fermions at L = 40 compared with total 
contribution, (c) Total contribution of the doublers to the anomaly integral for 0.01 < am < 0.1 
for naive Wilson fermions at L = 40 for 1 < r < 2. (d) Comparison of total doubler contribution 
and the physical fermion contribution for 0.01 < am < 0.1 with naive Wilson fermions at L = 40 
and r = ±1. 

L > 50. Nevertheless the infinte volume result for OStm fermion is closer to the continuum 
compared to 0(a) improved result. The small L behavior is different for the three types 
of fermions considered here. 

We have checked the sensitivity with respect to boundary condition of the volume 
dependence of the anomaly integral for the three actions in Fig. 3(d). For the 0(a) 
improved Wilson fermion, there is sensitivity to boundary condition only for very small L. 
On the other hand both for the naive and OStm Wilson fermion, the APBC improves the 
convergence of the anomaly integral with respect to L. 

In Figs. 3(a), 3(b), 3(c) and 3(d) we have seen some interplay between the finite volume 
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Figure 3. (a) Volume dependence of the anomaly integral for r = 1 and 0.01 < am < 0.1 for 
naive Wilson fermions with PBC. (b) Volume dependence of the anomaly integral for r = 1 and 
0.003 < am < 0.1 for 0(a) improved Wilson fermions with PBC. (c) Volume dependence of the 
anomaly integral for r = 1 and 0.003 < am < 0.1 for OStm Wilson fermions with PBC. (d) Volume 
dependence of the anomaly integral for am = 0.003 and 10 < L < 120 for naive Wilson, O(a) 
improved Wilson and OStm Wilson fermions with PBC and APBC. 

and the cutoff effects for naive Wilson fermion. The fact that this interplay is absent for 
0(a) improved Wilson fermion indicates that this is largely an 0(a) artifact. Finally Fig. 
4 shows that the convergence of the anomaly integral with respect to L is vastly improved 
by r averaging, since it removes the cutoff effects to a great extent. 



6 Pseudoscalar density term versus the anomaly term 

After the detailed analysis of the anomaly integral with respect to the am, r and L depen- 
dence, we revisit the flavor sing let axial Ward Takahashi identity, Eq. (2.1). To 0(g 2 ), for 
slowly varying background gauge fields the contribution from the PSD term is 
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Figure 4. Volume dependence of the anomaly integral for am = 0.003 and 10 < L < 130 for 
r = ±1 and r averaged for naive Wilson with PBC. 



2m (ipx^x) = 2 g 2 e M „ pA trace F flI/ (x)F pX {x) —^'^cos(p IJ ,a)cos(p 1 ,a)cos(p p a) 



x m 



cos(p x a)[m + Wo(p)] - 4 - sin 2 (p A o) {Go{p)) 

a J 

2 

-tt e^ P A trace F flu (x)F p x{x)I m (am,r, L) 



(6.1) 
(6.2) 



where 



I m (am,r,L) = cos(p jLt a)cos(p i/ a)cos(p p a) 



x m 



|cos(p A a)[' 



+ W (p)] - 4 - sin 2 (p A a) (S (p)) 3 ■ (6.3) 
a J 



First by using Reisz power counting theorem [20] we show that the PSD term to 0(g 2 ) 
is independent of m in the infinite volume continuum limit. In the limit L — > oo 



I m (am,r) 



4 N(m,p,a) 



D(m,p, a) 



(6.4) 
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where 



N(m,p,a) = cos(pua)cos(p u a)cos(p p a) x m cos(p\a)[m + Wq(p)] — 4 — sin 2 (p\a 



= cos(p„a)cos(pj,a)cos(p p a) x m cos(p\a)[m H — [1 — cos(ap„)l] — 4 - sin 2 (p^o) 

L a ^— ' a 

m 

D(m,p,a) = Go(p)~ 3 = I ^ ^sin 2 (ap A1 ) + (m + -^[l-cos(ap M )]) 2 j 

V M M / 

To find lattice degree of divergences (LDD) of N(m,p,a) and D(m,p,a) we scale p and a 
in the following way such that p = pa remains fixed, 

p — > Xp 
a — > a/A. 

Then, 

N(m,pX,a/X) = -[m^[l - cos(ap M )] - 4 sin 2 (p A a) A + O(A ) 

M 

r 1 3r 2 
D(m,pX,a/X) = U sin 2 (ap M )] 3 + -^£[1 - cos(a^)]] 6 + sin 2 (ap M )] 2 [^[l - cos(op^)]] 1 

L Hi Li Li 

MM MM 

+ l^[^ s in 2 (ap M )]E[l-cos(a^)]] 4 ] A 6 + 0(X 5 ) 

M M 

Then according to Reisz, LDD of N (degrpN) and LDD of D (degrpD) are given by, 

degrpN = 1 
degrpD = 6 

So, LDD of the integral in Eq. (6.4) is 

4 + 1- 6 = -K0 

Then according to the Reisz power counting theorem, continuum limit of I m (am,r) is 



POO 



d p lim cos(pfj,a)cos(p u a)cos(p f 



a) x m 



cos(p A a)[m + W (p)] - 4 - sin (p A a) ((?o(p)) c 



d 4 p 



7T" 



(p 2 + m 2 ) 3 2 



I m (am,r,L) = ~I m (am,r,L) . 



7T 



(6.5) 



lim I m (am,r,L) = —1 
a— >0, L— >oo 

and hence lim a _^o, L->oo Im(am, r, L) = — 1 is independent of m. 
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For the PSD term when am decreases one needs to go to larger value of L compared 
to the anomaly term to achieve the infinite volume limit. It is important to note that as 
long as am is nonzero, there is no infrared singularity. This is borne out by numerical 
calculations as shown in Fig. 5(b). Furthermore we see that the lattice result differs from 
the continuum result by exhibiting considerable am dependence (see Fig. 5(a)) which 
appears to be a manifestation of cutoff effects with naive Wilson fermion. Nevertheless 
we see that the two terms, namely, the anomaly term and the PSD term to 0{g 2 ) cancel 
each other independent of the value of am for large enough volumes. To check whether r 
averaging can remove the cutoff effects we study the am dependence of the r averaged PSD 
term. Fig. 5(c) shows that the PSD term is almost independent of am (for am <= .1) as 
in the continuum when r-averaging is performed. 

7 Conclusions 

In this work we have addressed four issues related to the one loop lattice calculation of 
chiral anomaly in the case of Wilson-like fermions. 

Some general remarks are in order. The weak coupling perturbative analysis may not 
be sensitive to some lattice artifacts that may be present in numerical simulations. On 
the other hand, some of the cutoff effects found in the perturbative analysis especially for 
the naive Wilson fermions could be removed by smoothing of gauge fields and smearing of 
hadronic operators in numerical simulations. 

We have used one loop lattice calculation of chiral anomaly to study quantitatively the 
effects of r averaging which has been proposed before to achieve better scaling behavior of 
observables with Wilson fermions. We find that r averaging of naive Wilson fermion has 
a slightly better approach to continuum chiral limit compared to OStm Wilson fermions. 
We show that in this case r averaged result is much better than O(a) and even than 0(a 2 ) 
improved Wilson fermion with tree level coefficients. It is the physical fermion contribution 
which is largely influenced by the r averaging as seen in Fig. 2(d). 

We have studied the doubler contributions as a function of r and lattice fermion mass 
am. We have verified that the doubler contribution decreases as r increases and it is not 
very sensitive to the lattice fermion mass. 

Next, we have studied the possible interplay between finite size and cutoff effects by 
investigating in detail naive, 0(a) improved and OStm Wilson fermion cases for a range 
of volumes and am. We find that the major contribution to the interplay between finite 
volume and cutoff effects is at the 0(a) level for naive Wilson fermion. 

Lastly we have studied the relative roles played by the PSD term and the anomaly 
term in the flavor singlet axial Ward Takahashi identity for naive Wilson fermion on the 
lattice at one loop level for slowly varying background gauge fields. First by using Reisz 
power counting theorem we have shown that the PSD term to 0(g 2 ) is independent of m 
and cancels the anomaly term in infinite volume continuum limit. The lattice result for 
PSD term differs from the continuum result by exhibiting considerable am dependence for 
r = +1 or r = —1 (see Figs. 5(a) and 5(c)) which appears to be a manifestation of cutoff 
effects with naive Wilson fermion. Nevertheless we find that the two terms, namely, the 
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Figure 5. (a) The am dependence of the anomaly term and the PSD term for L = 100 at 0(g 2 ). 
(b) The L dependence of the anomaly term, the PSD term and the sum for am ~ 0.1 0(g 2 ). (c) 
The PSD term for r = ±1 and the average, for 0.03 < am < 0.2 and L = 300 0{g 2 ). 



anomaly term and the PSD term to 0(g 2 ) cancel each other independent of the value of am 
for r = +1 or r = -1 for large enough volumes. To verify whether r averaging can remove 
the cutoff effects we investigate the am dependence of the r averaged PSD term. Fig. 5(c) 
shows that the PSD term is almost independent of am when r-averaging is performed. 

In the light of the results obtained above, it is tempting to try out r-averaging in lattice 
QCD simulations with naive Wilson fermions. However, problems related with unphysical 
UV fluctuations may potentially make the simulation difficult. Smearing of gauge links in 
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the molecular dynamics trajectories has been used successfully to damp these unphysical 
fluctuations [21, 22]. 
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